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Abstract 

We propose a new method to extract 1^4 from the ratio of the decay dis- 



tributions *g = g[B X d l + r]/*§[B -> X s l+l~ }. This ratio depends only 
on the KM ratio 1^4 with 15% theoretical uncertainties, if dilepton invari- 
ant mass-squared s is away from the peaks of the possible resonance states, 
J/i/j, ip', and etc. We also give a detailed analytical and numerical analysis on 
f [a-X g i+r](g = d, S )andg . 
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1 Introduction 



The determination of the elements of KM matrix is one of the most important issues 
in the quark flavor physics. Moreover, the element Vtd (or Vub) is especially important 
to the Standard Model description of CP-violation. If it were zero, there would be no 
CP- violation from the KM matrix elements (i.e. within the Standard Model), and we 
have to seek for other source of CP violation in Kl — > 7T7t. Here we study the ratio 
In the Standard Model with the unitarity of KM matrix, Vt s is approximated by 
—V c b, which is directly measured by semileptonic B decays. There are already several 
ways to determine V t d available in the literature: 

• \Vtd\ can be indirectly extracted through — Bd mixing. However, in Bd — Bd 
mixing the large uncertainty of the hadronic matrix elements prevents us to 
extract the element of KM with good accuracy. 

• A better extraction of l^l can be made if B s — B s is measured as well, because 

1 Vis 1 ' 

the ratio f 2 B D I f 2 B D can be much better determined. 

J B d B d l J B a B s 

• The determination of I Yr- 1 from the ratios of rates of several hadronic two-body 
B decays, such as T{B° -> K*°K°)/T(B -> <f>K°), T(B° -> K*°K*°)/T(B° -> 
<PK*°), T(B+ -> K*°K+)/V(B+ -> (f)K+), and T(B+ -> K*°K* + )/V(B+ -> 
cf)K* + ) has been also proposed in Ref. [|l|]. 

• V t d can be determined from K + — > 7i + uu decay within ~ 15% theoretical un- 
certainties with the branching fraction B ~ 10 -10 [0. 

Why do we discuss yet another method to determine Vt<p- The main reason we 
are interested in B physics is that this area is very likely to yield information about 
new physics beyond the Standard Model. We expect that new physics will influence 
experimentally measureable quantities in different ways. For example, most of us 
expect that AB = 2 transition is more sensitive to new physics than the decay rates. 
New physics may couple differently to K mesons compared to B mesons. Therefore, 
it is essential to determine the KM matrix elements in as many different methods as 
possible. 

In this paper, we propose another method to determine j^M precisely from the 
decay distributions ^[B — > X s l + l~ ] and ^jj[B — » where s is invariant mass- 

2 



squared of final lepton pair, . In the decays of B — > X q l + l~ (q = d,s), the short 
distance (SD) contribution comes from the top quark loop diagrams and the long 
distance (LD) contribution comes from the decay chains due to intermediate char- 
monium states. Therefore, the former (SD) amplitude is proportional to Vt q *Vtb, and 
the latter (LD) proportional to V cq *V c b. If the invariant mass-squared of is away 
from the peaks of the charmoniuum resonances {J/ip and ?/>'), the SD contribution 
is dominant, while on the peaks of the resonances the LD contribution is dominant, 
and therefore we expect that in the SU(3) symmetry limit the ratio 

f = f ( B^x d in-)/f(B^x,in-) 

as as as 



becomes 



1 dR(s) 

Y 2 ds 



pf + r, 2 (s^lGeV 2 ), 

1 (a-™,,/)' w 



where A is sin8 c , 8 C = Cabibbo angle, and for p, 7], see Eq. ([|). In the intermediate 
region, there is a characteristic interference beteen the LD contribution and the SD 
contribution, which requires the detailed study of the distributions. By focusing on 
the region, l(GeV 2 ) < s < rn 2 ,^, we can study |^| from the experimental ratio of 



the distributions 4^. 

as 

We stress the advantages to use the inclusive semileptonic decays: 

• If dilepton invariant mass-squared s is away from the peaks of the possible 
intermediate states, J/ip, if)', p, u>, and etc., the short distance contribution due 
to top quark loop is dominant. Therefore, we may extract very precisely the 
combination of |^| from the decay distributions within the range l(GeV 2 ) < 
s < mj/q, 2 without any theoretical uncertainties from the LD hadronic matrices, 
unknown KM matrix elements, and etc. 

• Inclusive decay distributions are theoretically well predicted by heavy quark 
mass expansion away from the phase space boundary. The leading order of the 
expansion agrees with the parton model result. Furthermore, non-perturbative 
1/nib 2 power corrections of QCD can be easily incorporated. 

This paper is organized as follows. In Section 2, we present the analytic formulae 
for ^[B X s l+l~ ] and ^[B -> X d l + l~\. We also show in detail the decay distri- 



bution ^r[B — ► X s l + l~ ] including the uncertainties coming from top quark mass m t 
and the scale of renormalization group /i. In Section 3, after brief discussion on the 
limit of KM elements coming from Bd — Bd mixing and B — > X u lT>, we study the ratio 
of the decay distributions, and show how we may extract \yj L \- 

2 Effective Hamiltonian for b — > ql + l~ and their 
differential decay rates. 

In this Section, the differential decay rates for B — > X q l + l~ (q = d, s) are shown 
including both LD and SD effects as well as Kq C d 2 l m b 2 power corrections. (See || 
for b -> dl + l~, @ for 6 -> sZ+Z", and §, § for B -> X s /+Z~ including the l/m b 2 
power corrections.) We also show in detail how the differential decay rate varies by 
changing the input parameters m t , and the renormalization scale /x. In Table 1, we 
summarize all the values of the input parameters used in our numerical calculations 
of decay rates. We use the central values for those input parameters, unless otherwise 
specified. 

The effective Hamiltonian for b — > ql + l~ (q = d, s) is given as 



■ 10 

EGA 

i=i 



AC 

+ -^V: q V ub [C l {0^-O l )+C 2 {0^-0 2 )\, (2) 
where V%j are the KM matrix elements. The operators are given as 

Ox = {q~La r yfibLa){cLp r y' J 'C L p), 

3 = {q La l^b La ) Y Wlp^Ilp), 

q'=u,d,s,c,b 
q'=u,d,s,c,b 

5 = (q La l^ La ) Y (?VT%/j)> 

q' =u,d,s,c,b 

Oq = (q La lvb L( ]) Y (eWAftJ. 

q'=u,d,s,c,b 

g 

7 = Y^qaa^irribR + mqL^aF^, 
Os = ^q a T^a, u (m b R + m q L)b p G a ^, 
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0\ = (qLalnbLodiuLpl^ULp), 

e 2 
e 2 

Oio = j^qa^LbJ-f^l, (3) 

where L and R denote chiral projections, L(R) = 1/2(1 =p 75). We use the Wilson 
coefficients given in the literature (see, for example, 0). 

With the effective Hamiltonian in Eq. (^), the matrix element for the decays 
b — > ql + l~ (q = d, s) can be written as 

V2n 



M{b^ q m-) = -^v*v tb 



CSf-Cio) (q^Lb) (T-fLl) 



-2Cf (qi <v ^(m q L + m bJ R) foj (V Z 



(4) 



where Cf 1 ^ is given by 



C 9 f (S) = C 9 |l + ^w(S)} + + W(5) , (5) 

where q v is the four momentum of s = q 2 , and s = s/m 2 . The function YsD q {s) 
is the one-loop matrix element of O9, and 1ld 9 (s) is the LD contributions due to the 
vector mesons mesons J/ip, ip' and higher resonances. The function uj(s) represents 
the O(ot s ) correction from the one-gluon exchange in the matrix element of Og, and 
is given in our Appendix. The two functions Yso q and Yld 1 are written as 

Y SD \s) = g(m c ,s)(3C 1 + C 2 + 3C 3 + C 4 + 3C 5 + C e ) 
~g(l,s) (4C 3 + 4C 4 + 3C 5 + C 6 ) 

~g{0,S) (C 3 + 3C 4 ) 

+l(3C 3 + C 4 + 3C 5 + C 6 ) 

+ C 2 )(g(0,s)-g(m c ,s)) , (6) 







v t *v tb v * 


3 

~ K 

cr 


( v; q v cb 


{ v t *v tb 


X 


E 



Y LD \s) = 4« {-ISLZC®- ^|^(3C3 + C 4 + 3 C 5 + C 6 

v t * q v tb 

irrjK -> /+/-) My. 
r mi,) <; s ) .U r 2 - sm b 2 - i.\l v V v . ' 



(7) 



where the function g(m c = s), g(l, s), and g(0, s) represent c quark, b quark and 
u,d,s quark loop contributions, respectively. The two functions g(z,s) and g(0,s) 
are given in our Appendix. In Eq. (0), C (0) = 3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 , 
and the first two terms are dominant, as can be seen from the Table 2. 

It is convenient to write the relevant combinations of KM in terms of the Wolfen- 
stein parametrization. In the following, in addition to A and A, we choose (1 — p) 2 + r] 2 
and 0i = arg(-^r^fr) as independent variables, then we have 

V*V tb ~ X 2 A , 

§^ - 0(A 2 ), 
V td V tb ~ A 3 A(1 - p + ft,) , 

i 

exp(-i^i) 



1 - p) 2 + ?7 2 

KdKfe _ p-iy 
V t * d V tb l-p + iri' 



-1+ .S^ . (8) 



;i-p) 2 +77 2 

In the SD contribution of I? — * , the w-quark loop contribution is neglected 

due to the smallness of the combination V* s V ub compared with V* s V cb ~ — V t * s V tb , while 
in B — > X d l + l~ the term which is proportional to V* d V ub is maintained. In the LD 
contribution, there is a contribution coming from the gluonic penguin amplitude which 
is proportional to v" 6 (3^3 + C4 + 3C$ + Cq). This contribution is neglected because 
of the smallness of the Wilson coefficient (3C3 + C4 + 3C5 + Cq) compared with C^°\ 
We adopt k = 2.3 |§ to reproduce the rate of the decay chain B — > X s J/tp — > X s l + l~ . 
Note that the data determines only the combination, kC' ' = 0.88. 

By combining both SD and LD contributions as well as non-perturbative l/m b 2 
power corrections, the differential decay rate for B — > X q l + l~ (q = d, s) becomes 



■'2 

-u{s, m 9 )((l - m 2 ) 2 + s(l + m 2 q ) - 2s 2 ) 



dB 2 \v; q v tb 

\V cb \ 2 

H — (1 — 4m 2 . + 6rh q — 4m q + rh q — s + rh 2 s + m^s — m q s 
3 



6 



2s3 



—3s 2 — 2m 2 s 2 — 3m 4 s 2 + 5s 3 + 5m^s 



2s 4 



Ai 



(l — 8m 2 + 18m 4 — 16rh q + 5m q — s — 3m 2 s + 9m 4 s — 5m^s — 15s 

A 2 



-18m 2 q s 2 - 15mp 2 + 25s 3 + 25m 2 s 3 - 10s 4 



u[s,m q ) 



'10 



+ 



s,m g ) (2(1 + m 2 .)(l - m 



2\2 
1> 



1 + Um z q + rff q )s - (1 + rh z q )s' 



4 



+-(2 - 6m 2 + 4m 4 + 4mJ - 6m J + 2m J° 
—5s — 12m 2 s + 34m 4 s — 12m^s — 5m^s + 3s^ 



+29m^s^ + 29m*s 2 + 3m°s 2 + s c 



10m ? s + m q s — s — m q s 



,4-3 



2;4\ 



Ai 



u{s,m q ) 



+4 (-6 + 2m 2 + 20m 4 - 12mJ - Um 8 q + lOmf + 3s + 16m 2 s + 62mjs 
-56mjs - 25mjs + 3s 2 + 73mjs 2 + lOlmJs 2 + 15mjs 2 + 5s 3 - 26mjs 3 + 



5m g s — 5s — 5m g s 



A 2 



u{s } m q ) 



eff 12 



yu^s, m g j 



m 



2\2 



-(1 + m g )sj + 4(1 - 2m g + m q - s - m q s) u(s, rh q ) Ai 

+4 (-5 + 30m 4 - 40mJ + 15m s q - s + 2\m 2 q s + 25m 4 s - 45mjs + 13s 2 



+22m 2 s 2 + 45m 4 s 2 - 7s 3 - 15m 2 s 3 



We explain some aspect of the expression, 



A, 



ite(cf) cf\ 



(9) 



The branching ratio is normalized by B s i of decays B — > (X c , X u )lvi. We sep- 

ly* |2 

arate a combination of the KM factor ^ ^ due to top-quark loop from the 
normalization factor £>o- The normalization constant Bq is 

3 a 2 1 



B n = B 



si: 



(10) 



16tt 2 /(m c )«(m c ) ' 

where /(m c ) is a phase space factor, and /«(m c ) accounts for both the 0(a s ) 
QCD correction to the semi-leptonic decay width and the leading order (l/m?,) 2 
power correction, They are given in our Appendix explicitly. 



The Wilson coefficients depend on the top quark mass m t , the renormalization 
scale fi, and Aq C d- Their dependences are studied in Tables 2, 3 and 4. As 
the value of \i (Aqcd) decreases (increases) , |CV| and C$ are getting larger, 
while C10 is independent on //. As m t increases, |CV| , Cg and |C 10 | all become 
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larger. The matrix elements also depend on the renormalization scale /x, and 
their dependence is partially cancelled by the given dependence of the Wilson 
coefficients. The overall dependence can be studied with the differential rate. 
In Figs. [I] and 0, we show the dependence of the differential rate on m t and 
/i. As the value of m t increases from 166 (GeV) to 184 (GeV), the differential 
rate increases about 20% at s ~ 5 (GeV 2 ), as shown in Fig. [I]. At around the 
same region of s, by changing \i from 2.5 (GeV) to 10 (GeV) the differential rate 
increases about 20%, as shown in Fig. |2|. This observation is consistent with 
the result of Ref . [7j , in which only the SD contribution has been analysed. 

The differential decay rate @ is not a simple parton model result. It contains 
non-perturbative \jvn 2 power corrections, which are denoted in Eq. (|J) by 
the terms proportional to Ai and X 2 - The parameters Ai and A2 are related to 
the matrix elements of the higher derivative operators of heavy quark effective 
theory §, @; 



B 



h(iD) 2 h 



B 



h —a^ G, v h 
2 M 



2M B Ai = 2M B m b 2 X 1 
6M B X 2 = 6M B m b 2 X 2 



(11) 



where B denotes the pseudoscalar B meson, is the covariant derivative, and 
is the QCD field strength tensor. 



3 The extraction of 



Yu 
v ts 



Presently we can constrain the value of (1 — p) 2 + rj 2 from B& — B& mixing, while from 
B — > X U W we obtain some limit on p 2 + rj 2 . Using those given limits, we can show 
how the ratio w^r depends on input KM values, (1 — p) 2 + rj 2 and 0i (or (3). 

As is well known, |V^V^| are written with AM^ of Bd — B^ mixing. It reads as 



\V t tVu\ 



AX' 



1-p) 2 
6tt 2 



1 l/2 



G 2 F % CD M^M B \ 



AM, 1 ' 2 



Blt 2 f B \E{x t )\V* 



rat 



(12) 

2- It 



where E(x t ) is the Inami-Lim function ||TD[ of the box diagram, and x t = Ti 
has the value |S(xt)| = 2.58 for m t = 175 (GeV) and it varies from 2.38 to 2.78 as 



S 



m t varies from 166 (GeV) to 184 (GeV). We use the following values; rj QCD = 0.55, 
m t = 175 (GeV), and the experimental constraints 

AM d = 0.474 ±0.03lG»s~ 1 ) , 



cb 



Xjp 2 + V 2 = 0-08 ±0.02 . (13) 



Then, those constraints (13) are translated into the limits of (1— p) 2 +r) 2 and a/ p 2 + rj 2 



as 

(l-pf + ri 2 = 0.85(1 ±0.15) {0.2/(^/ 2 / fld (GeV))} 2 

= [0.59 ±0.09, 1.3 ±0.20] , 

y + r] 2 = o.36 ± 0.09 . (14) 

Here we used B)j 2 f B = 0.2 ± 0.04(GeV) to get the range of (1 - p) 2 + r] 2 « [0.5, 1.5]. 
We also used AX 2 = 0.041(±0.003), and A = 0.2205. In Fig. |, we show the present 
limit in the plane of (0i, (1 — p) 2 + i] 2 ). The horizontal axis corresponds to 4>i (or (3), 
and the unit is degree. The vertical axis corresponds to (1 — p) 2 + rf. The thin dashed 
line is obtained from the central values of \V u b\ and the thin solid line is obtained from 
\Vtd\- The central value corresponds to (0i, (1 — p) 2 + rj 2 ) « (20°, 0.85). 

Now let us consider the ratio A- dR } s ' . In the SU(3) limit, we expect that this 
ratio approaches to the input value of (1 — p) 2 + rj 2 for the dileptonic invariant mass- 
squared s far below the peak of charmonium resonances. If s is on the peak of 
resonances, the ratio becomes 1. In Fig. £|, we show the ratio dR !> s ' for two sets of 
the assumed input values of (0i, (1 — p) 2 + r/ 2 ), one of whose (1 — p) 2 + rj 2 is 0.59, 
and the other is 1.33. They correspond to their small (or large) values allowed from 
the present experimental result of Bd — mixing. The solid curve corresponds to 
(0i, (l-p) 2 + ?7 2 ) = (20°, 0.59), and the dot-dashed curve corresponds to (20°, 1.33). 
The ratios for the CP conjugate process B — > X q l + l~ are also shown, denoted by the 
long-dashed curve for (1 — p) 2 +rf = 1.33, and by the dashed curve for (1 — p) 2 +rf = 
0.59. They are obtained by reversing the sign of (pi in the corresponding B — * X q l + l~ 
process; i.e., 4>i — > —<j>i- They are labeled as (—20°, 1.33) and (—20°, 0.59) in Figure 
I 

To summarize the numerical results of Fig. ^ and Section 3: 

• The predicted ratio at low invariant mass region s ~ 1 (GeV 2 ) is very near to 
the our assumed input value of (1 — p) 2 + rf, while on the peak of the resonances 
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J/ip, if}', this ratio becomes almost 1, as we expected, i.e. 
1 dR(s) 



A 2 ds 



1 (s ~ mj/^ ,m^/ , ..) , 

(1 — p) 2 + rf (s being away from rrij/^ 2 , ..) . (15) 



In the intermediate region, there is a characteristic interference between the 
LD contribution and the SD contribution, which can be only derived from the 
detailed expression of the distributions, Eq. @. 

The value of the ratio does not depend much on whether the decaying particles 
are B or B at any invariant mass-squared region. 

This ratio changes only a few % when we change the input parameters, m t and 
/i, within the range shown in Table 1. The dependences on mt and \x are almost 
cancelled away in the ratio. Therefore, the uncertainties in this ratio due to 
the input parameters are much smaller than the uncertainties in the differential 
decay rate itself, as shown in Tables 2, 3 and 4 as well as in Figures [l] and [|. 

In Fig. [| the range between the dot-dashed curve and the solid curve corre- 
sponds to the value of the ratio tj^^ allowed from the present experimental 
result of Bd — Bd mixing. Future experimental measurements on the ratio of 
the branching fractions, B(B — > XJ, + l~) and BiB — ► XJ + l~), can give much 
better alternative for determination of |^| without any hadronic uncertainties, 
limited only by experimental statistics. 

If 10 9 BB pairs are produced, the expected number of the events for B — > 
XJ + l~ in the range of 4m M 2 < s < 6 GeV (m M = muon mass) is about 100 
for (1 - pf + 7] 2 = 0.59, and is about 220 for (1 - p) 2 + rf = 1.33. Therefore, 
the statistical accuracy of (1 — p) 2 + rf determined from this method is about 
7% ~ 10% with the expected production of 10 9 BB pairs. 

We have assumed in our numerical analysis the flavor SU(3) symmetry with 
m d = m s- We estimated the corrections due to SU(3) breaking by varying m d 
from 0.01 GeV up to m s = 0.2 GeV. And we find the ratio decreases within 
0.2 ~ 0.3% for the range 1 < s < 9 (GeV 2 ). Therefore, we conclude the SU(3) 
breaking effect does not affect the extraction of at all. 

° 1 Vts 1 
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Parameter 


Value 


m w 


80.26 (GeV) 


m z 


91.19 (GeV) 


• 2 n 

sin Uw 


0.2325 


m s 


O.Z [Kje\ ) 


m d 


U.U1 ^ Vjre V J 


m c 


1.4 (GeV) 


m b 


4.8 (GeV) 


m t 


175 ± 9 (GeV) 


n 




a(5) 
^QCD 


0.214lg;8g (GeV) 


a QED 


129 


a s (m z ) 


0.117 ±0.005 


Bsi 


(10.4 ±0.4) % 
-0.20 (GeV 2 ) 
±0.12 (GeV 2 ) 


Ai 


A 2 



Table 1: Values of the input parameters used in the numerical calculations of the 
decay rates. Unless, otherwise specified, we use the central values. 



n 


Ci 


c 2 


c 3 


c, 


c 5 


c 6 


cf 


qNDR 


c(°) 


5 


-0.2404 


1.1031 


0.0107 


-0.0249 


0.0072 


-0.0302 


-0.3110 


4.1530 


0.3805 


10 


-0.1606 


1.0642 


0.0068 


-0.0170 


0.0051 


-0.0194 


-0.2768 


3.7551 


0.5816 


2.5 


-0.3472 


1.1614 


0.0163 


-0.0348 


0.0096 


-0.0462 


-0.3525 


4.4128 


0.1163 



Table 2: fj, (in GeV) dependence of the Wilson coefficients used in the numerical 
calculations. The values of Aqcd and m t are fixed at their central values. 
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Aqcd 5 


c x 


c 2 


c 3 


c A 


c 5 


c 6 


C 7 etf 


(jNDR 


c(°) 


0.214 


-0.2404 


1.1031 


0.0107 


-0.0249 


0.0072 


-0.0302 


-0.3110 


4.1530 


0.3805 


0.280 


-0.2579 


1.1122 


0.0116 


-0.0265 


0.0076 


-0.0327 


-0.3181 


4.2137 


0.3369 


0.160 


-0.2242 


1.0949 


0.0099 


-0.0233 


0.0068 


-0.0279 


-0.3043 


4.0891 


0.4212 



Table 3: Aq C d 5 dependence of the Wilson coefficients used in the numerical calcu- 
lations. The values of m t and fi are fixed at their central values. 



m t 




(jNDR 


Cio 


175 


-0.311 


4.153 


-4.546 


166 


-0.3066 


4.0796 


-4.1877 


184 


-0.3150 


4.2238 


-4.9156 



Table 4: m t dependence of the Wilson coefficients used in the numerical calculations. 
The values of Aqcd 5 and \i are fixed at their central values. 



Appendix 



A Functions g(z,s), g(0,s) and uj(s 



The functions g(z, s) and g(Q, s) are given as 

,m 6 , 8. 8 4 



9\z, s) 



9 ln(^)--ln, + - + -y--(2 + y )V^[ 



27 9 l 



x 



8(1 - y) (In ] + ^ } V - m) + 9(y - 1)2 arctan ' 



with y = 4z 2 /s, and 



9(0, s) 



8 8 ,m b 4 4. 

In — Ins H — ?7r 

27 9 v /i 9 9 



(16) 



(17) 



The function u)(s) represents the 0(a s ) correction from the one-gluon exchange 
in the matrix element of Og 



OJIS 



jvr 2 - iu*{s) - \ In51n(l - s) - J*,^ ln(l 



2s(l + s)(l-2s) 
3(1 - s) 2 (l + 2s) 



In s + 



3(1 + 2s) 
5 + 9s - 6s 2 
6(1 - s)(l + 2s) ' 



;is) 
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B Functions f(m c ) and /c(m c ) 



The phase space function for T(B — > XJu) in the lowest order (i.e., parton model) 



f(m c ) = 1 — 8 ml + 8 ml — — 24 m\ mm c 



(19) 



And K(m c ) accounts for both the 0(a s ) QCD correction to the semi-leptonic decay 
width and the leading order (1/m;,) 2 power correction, 



K(m r 



2a s (m b ) A h(m 

l = g(m c ) + 



3tt 



2m\ 



(20) 



The function g(m c ) is given [12], [13 

g{rh c ) 



as 



7r 2_H)(l_ 7 y lc )2+3 

4 ' 2 



(21) 



and finally the function h(m c ) is given || as 
A 2 



h(m c ) = Ai + 



/(m c 



9 + 24m^ - 72m* + 72m° - 15m° - 72m* lnm c 



(22) 



where Ai (or — /i^) and A 2 (or /i^) denote the matrix element of the higher derivative 
operators of heavy quark effective theory, as defined in Eq. ([□]) . 
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10.0 



mt=175 




0.0 5.0 10.0 15.0 20.0 25.0 

S ( GeV 2 ) 



Figure 1: Dilepton invariant mass spectrum, ^jf (-B — > X s e + e~) x 10 7 . The unit of 
s is (GeV 2 ). The solid line corresponds to m t = 175 (GeV). The short-dashed line 
corresponds to m t = 175 + 9. The long-dashed line corresponds to m t = 175 — 9. 
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15.0 




-11.0 - 
-13.0 - 



0.0 5.0 10.0 15.0 20.0 

S (GeV 2 ) 

Figure 2: The % variation of dilepton invariant mass spectrum defined as ( — 
^)/^; B = B(n = 5(GeV)). The solid line corresponds to ^ = 10 and the dashed 
line corresponds to /i = 2.5. 
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2.0 - 



B d B d bar Mixing 

V„ h 



1.6 



+ 12 - 



0.8 - 



0.4 



0.0 



0.0 5.0 



10.0 



15.0 20.0 
^ (degree) 



25.0 30.0 



Figure 3: The limit on (1 — p) 2 + if and 0i (or (3). The horizontal axis corresponds 
to 0i, and the unit is degree. The vertical axis corresponds to (1 — p) 2 + if. The 
thin dashed line and the thin solid line are obtained from the central values of 1 
and \Vtd\ respectively. The thick solid lines are obtained from the allowed range of 
\Vtd\ from BdBd mixing. The thick dashed lines are obtained from the allowed range 
of IK&I- 
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Figure 4: The ratio versus s (GeV 2 ) with two different input values of (1 — 

p) 2 + r\ 2 . The solid curve correspond to (fa, (1 — p) 2 + if) = (20°, 0.59), and the 
dot-dashed curve corresponds to (20°, 1.33). The ratios for CP conjugate process 
B — > X q l + l~ are denoted by the dashed curve for (1 — p) 2 + rf = 1.33, and by the 
long-dashed curve for (1 — p) 2 + rj 2 — 0.59. They are obtained by reversing the sign 
of 0i in the corresponding B — > X q l + l~ process; i.e., <pi — > —fa. 
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